Kondo screening in unconventional superconductors: 
The role of anomalous propagators 
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The Kondo effect in superconductors is frequently investigated using the local quasiparticle density 
of states as sole bath characteristics, i.e., the presence of anomalous propagators is ignored. Here we 
point out that this treatment is exact for a number of situations, including point-like impurities in d- 
wave superconductors. We comment on recent investigations [M. Matsumoto and M. Koga, J. Phys. 
Soc. Jpn. 70, 2860 (2001) and Phys. Rev. B 65, 024508 (2002)] which reached different conclusions: 
while their numerical results are likely correct, their interpretation in terms of two-channel Kondo 
physics and an "orbital effect of Cooper pairs" is incorrect. 

PACS numbers: 75.20.Hr,74.70.-b 



The physics of quantum impurity moments in super- 
conductors, associated with the Kondo effect, has been 
subject of numerous investigations in recent years. Di- 
verse theoretical techniques have been employed to study 
Kondo or Anderson models in a superconducting envi- 
ronment [1-9] . Most of these studies effectively neglect 
the presence of superconducting (sc) fluctuations in the 
host, i.e., they use the local fermionic density of states 
(DOS) as the only input quantity characterizing the en- 
vironment of the Kondo impurity. For s-wave supercon- 
ductors it has been shown [1] that this approximation 
is not justified: here the properties of a Kondo impu- 
rity are different from the ones of an impurity embedded 
in a non-superconducting system with the same DOS. 
More precisely, the superconducting bath turns out to 
be equivalent to a non-superconducting bath with addi- 
tional particle-hole (p-h) asymmetry [1] . As a result, for a 
p-h symmetric conduction band a screened singlet state 
is realized at large Kondo coupling, in contrast to the 
non-sc case with a hard gap in the local DOS [11]. This 
difference can be understood as caused by the anomalous 
bath propagators. 

In this Brief Report, we address the role of anoma- 
lous propagators in unconventional superconductors. We 
argue below that neglecting sc propagators is exact in 
many, potentially experimentally relevant, cases, e.g., for 
point-like impurities in d-wave and unitary p-wave su- 
perconductors. In these situations, the dynamics of the 
impurity degrees of freedom can be calculated using the 
normal bath propagators only. We also comment on re- 
cent papers by Matsumoto and Koga who argued in favor 
of a non-trivial "orbital effect of Cooper pairs" for point- 
like impurities in both p + ip and d + id superconductors 
[7, 8]. While we believe that their numerical results are 
correct, we point out that such an orbital effect does not 
exist: in their situation the impurity properties are ex- 
clusively determined by the local DOS. 

We start from the action of an Anderson impurity in 



a general interacting host. 



S = j^^2^2c kcr (iuj n ){-iuj n + e k ]c ka (iuj n ) (1) 



uJn ko 



f 

S- mt (c,c) +<Sioc(/a) + / dr yjCVj/gCjg + c.c.) 



{i}<r 



where (3 = l/T is the inverse temperature, c are the 
conduction electrons with dispersion e k on a regular lat- 
tice with Af sites, <Si n t are the interactions within the 
conduction band, and <5>i oc describes the / electron im- 
purity orbital with on-site energy and repulsion. The 
sum runs over a set of lattice sites in the vicin- 

ity of the impurity, and Vi is the hybridization matrix 
element; for a point-like impurity only a single Vi is non- 
zero. We can define a linear combination cq of conduc- 
tion electron operators to which the impurity couples: 
Vc a = Y,{i} v i c i<? witn [ C 0;C,! ) ]+ = 1. After decoupling 
of <S; nt and a suitable saddle-point approximation of BCS 
type, all conduction electrons except cq can be integrated 
out. (After the BCS approximation, the remaining in- 
tegral is Gaussian and can be performed exactly.) To 
simplify notation we will restrict ourselves to BCS sin- 
glet pairing; the arguments apply similarly to unitary 
triplet states, whereas non-unitary triplet states require 
an additional coupling between the impurity and the 
condensate spin moment. Introducing a Nambu spinor 
= (coti c o|) we obtain an action of the form: 

S = -— ty (iuj n ) G^jiLUn) yr (iu) n ) 



+ Sioc(U) + / dr(V f a c 0a + c.c.) 



(2) 



Here, Go(iuj n ) is the local conduction electron Green's 
function at the impurity location. In principle, a treat- 
ment of iSint beyond mean-field also generates a retarded 
sclf-interation for the cq. As we are mainly interested 
in a BCS-type host we will neglect this here. Then, the 
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FIG. 1: Spatially extended impurities for which the anoma- 
lous propagator in Go drops out, i.e., ~^2 k hj.Ak/Dk = 0, for 
a square-lattice d-wave superconductor. The dot represents 
the impurity orbital, the dashed lines show the non-zero hy- 
bridization paths (in the notation of an Anderson model). 
The situations b) and c) correspond to the s-wave and d-wave 
linear combinations, i.e., screening channels, of the four host 
sites, which occur in a Kondo model for a spatially extended 
impurity, as discussed in Ref. 5. 

properties of the bath are completely contained in Go- 
Explicitly we have 

where Ak is the complex gap function, and the function 
hk contains the geometry of the impurity coupling to the 
host: 

h k = Y,e ikRi V i /V. (4) 
{<} 

For a point-like impurity hk = 1, and the anoma- 
lous (off-diagonal) part of Go is given by Gq(^„) = 
J2 k Ak/Dk (with D k = u>n + el + |A fe | 2 ), which vanishes 
for unconventional superconductors with inversion sym- 
metry and Cooper pair angular momentum I > 0, i.e., 
p-wave, d-wave or higher symmetries. (This also applies 
to p x + ip y or d x 2_ y 2 + id xy pairing states, but not neces- 
sarily to d+is states.) Then, only the diagonal part of Go 
enters the impurity action (2), which is completely deter- 
mined by the local DOS. For a spatially extended impu- 
rity hybridized with more than one site the anomalous 
piece of Go is given by ^2 k h\AkjDk which still vanishes 
for some important situations: Consider a d-wave super- 
conductor, Ak = Ao(cosfc x — cosky), and an impurity 
hybridized e.g. with four sites as in Figs. la,b,c. In all 
cases the average over h\Ak/Dk vanishes for symmetry 
reasons. 

At this point a brief comment on the experimental sit- 
uation is in order. Point-like impurities can be realized in 
layered superconductors with out-of-plane impurity mo- 
ments coupling to a single conduction electron orbital 
only. In contrast, in-plane impurities will typically cou- 
ple to a number of sites, see Fig. 1. The signs of the 
various hybridization matrix elements depend on the in- 
volved orbitals of both impurity and host atoms. We note 
that the situations in Figs. lb,c were used in Refs. 4, 5 to 
model the magnetic moment induced by a Zn impurity 
in a high-temperature superconductor. (There, a Kondo 
model for a spatially extended impurity was employed, 
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FIG. 2: Numerical results for the local susceptibility, ob- 
tained from a single-band NRG for a (semiconducting) bath 
with the hard-gap DOS of a d + id (or p + ip) superconduc- 
tor. The host bandwidth is unity, the NRG parameters are 
A = 3 and N s — 1000. a) Txioc for different gap values 
A = 10~ 2 , 10~ 3 , 10~ 4 , 10~ 5 (from top to bottom) where the 
Kondo temperature Tk of the model at A = is 3 x 10~ 3 . 
b) The value of the local Curie constant C, normalized to its 
free-spin value 1/4, as a function of the ratio Tk/A. Differ- 
ent symbols correspond to different values of A. (Deviations 
from a single universal curve are primarily due to NRG dis- 
cretization effects.) The results can be compared with Fig. 4 
of Ref. 8. 

leading to multiple screening channels. In addition, the 
situation for Zn is complicated by the fact that Zn, having 
a filled d shell, acts as a vacancy, but induces a magnetic 
moment in its vicinity - for details see Ref. 5.) 

Returning to the models discussion - what about 
the behavior of a spatially extended Anderson impurity, 
where Gg does not vanish? Formally, we still have a 
single-channel model (2), but now in the presence of 
both a normal and an anomalous bath - this is sim- 
ilar to a point-like impurity in a s-wave superconduc- 
tor. As explained in Refs. 1, 12 the main effect of the 
anomalous bath can be understood as a transverse charge 
pseudospin field which induces an additional particle-hole 
asymmetry. 

So far we have discussed Anderson impurity models. 
For a single-site impurity the above discussion applies 
identically to a Kondo model (via Schricffcr- Wolff trans- 
formation). Multichannel physics arises only in a spa- 
tially extended Kondo impurity model (independent of 
anomalous propagators), see Refs. 5, 13. 

Now we comment on recent papers by Matsumoto and 
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Koga [7, 8] about a non-trivial "orbital effect of Cooper 
pairs" on the dynamics of point-like impurities in p + ip 
and d + id superconductors. They argue that, although 
the impurity only couples to the s-wave channel of the 
conduction band, the I ^ Cooper pairs mediate an 
indirect coupling to higher angular momentum channels, 
leading eventually to a multichannel Kondo problem. For 
the cases of p + ip and d + id pairing, where the problem 
reduces to two angular-momentum channels, they pro- 
vided a numerical solution utilizing a two-band numeri- 
cal renormalization group (NRG) method. Their findings 
show the absence of screening even for large Kondo cou- 
pling, in contrast to the s-wave situation. We believe that 
their mapping and results are likely to be correct, but al- 
low for a much simpler interpretation. Specifically: (i) 
Our above arguments, which are rigorous for a point- like 
impurity in a c?-wave or unitary p-wave BCS supercon- 
ductor, show that only the local DOS enters the impu- 
rity problem, (ii) The two-band model of Ref. 8 [their 
Eq. (2.21)] is in fact equivalent to a one-band model of 
the form (2) (in that it yields the same partition func- 
tion and the same observables related to the impurity 
degrees of freedom). A first indication is that the impu- 
rity term in Eq. (2.21) of Ref. 8 mixes the two bands, i.e., 
contains channel-flip terms (in contrast to true two-band 
models). Formally, a one-band model can be obtained 
by integrating out the dk,i=i band in Eq. (2.4) of Ref. 8 
exactly - this will generate a normal self-energy for the 
fl/U=o fermions and induce the quasiparticle gap [10]. (iii) 
The local DOS of a p + ip or d + id superconductor dis- 
plays a hard gap, thus we expect the qualitative impurity 
properties to be those of a hard-gap system [11]. In par- 
ticular, for a p-h symmetric band there is no screening 
even at large Kondo coupling, consistent with the results 
of Ref. 11. However, in the presence of p-h asymmetry 
there will be a first-order transition to a screened phase 
as the Kondo coupling is increased [11]. This transition 
was not found in Refs. 7, 8, instead the authors claim 
that their findings persist in the presence of p-h asym- 
metry (but numerical results are not shown) . We believe 
this claim to be incorrect, (iv) We have performed one- 
band NRG calculations using the local quasiparticle DOS 
of a d + id (or p + ip) superconductor, and have essen- 
tially reproduced the results for the local susceptibility 
of Ref. 8, see Fig. 2 [14]. Taken together, the physics of 
point-like impurities described in Refs. 7, 8 is completely 
contained in a Kondo model where the impurity is cou- 
pled to a single band with a (fully gapped) local DOS of 
a non-superconducting host [10], and the "orbital effect 
of the Cooper pairs" does not exist. (This is different for 



spatially extended Kondo impurities where true multi- 
band effects obtain [9].) The discussion in Sec. Ill B of 
Ref. 8 is formally correct, but misses that the mechanism 
leading to the doublet ground state of a strongly coupled 
Kondo impurity is identical for any p-h symmetric hard- 
gap system [11]. 

In summary, we have argued that in unconventional 
BCS superconductors the properties of a Kondo impurity 
are not influenced by the presence of anomalous propa- 
gators, provided that J^ fe h^Ak/Dk = (4) - this is e.g. 
the case for point-like impurities and a vanishing local su- 
perconducting order parameter. Under these conditions, 
using only the local density of states as bath input quan- 
tity for impurity calculations, as done in Refs. 4, 5, is 
exact. 
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